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PREFACE 


The  prototype  tests  described  in  this  report  were  conducted  during 
October  1980  by  the  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES)  under  the  sponsorship  of  the  U.  S.  Army  Engineer  District,  Seattle. 

Acknowledgment  is  made  to  individuals  of  the  Seattle  District  for 
their  assistance  in  the  investigation.  Mr.  E.  D.  Hart,  Chief  of  the 
Prototype  Evaluation  Branch,  was  test  coordinator  for  WES.  This  report 
was  prepared  by  Mr.  Hart  with  assistance  from  Mr.  E.  B.  Pickett, 

Mr.  J.  E.  Hall,  and  Dr.  F.  M.  Neilson,  under  the  general  supervision  of 
Mr.  M.  B.  Boyd,  Chief  of  the  Hydraulic  Analysis  Division,  and  Mr.  H.  B. 
Simmons,  Chief  of  the  Hydraulics  Laboratory,  WES. 

Commanders  and  Directors  of  WES  during  the  investigation  and  the 
preparation  and  publication  of  this  report  were  COL  Nelson  P.  Conover, 

CE,  and  COL  Tilford  C.  Creel,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con¬ 
verted  to  metric  (SI)  units  as  follows: 


_ Multiply _ 

acre-feet 

cubic  feet  per  second 
Fahrenheit  degrees 
feet 

feet  per  second 

feet  per  second  per 
second 

inches 

inches  per  second 

miles  (U.  S.  statute) 

pounds  (force)  per 
square  foot 

pounds  (force)  per 
square  inch 

pounds  (mass)  per 
cubic  foot 

slugs  (mass)  per 
cubic  foot 

square  feet 


_ & _ 

1233.482 

0.02831685 

5/9 

0.3048 

0.3048 

0.3048 

25.4 

25.4 

1.609344 

47.88026 

6894.757 

16.01846 

515.3788 

0.09290304 


_ To  Obtain _ 

cubic  metres 
cubic  metres  per  second 
Celsius  degrees  or  Kelvins* 
metres 

metres  per  second 

metres  per  second  per 
second 

millimetres 

millimetres  per  second 

kilometres 

pascals 

pascals 

kilograms  per  cubic  foot 
kilograms  per  cubic  metre 
square  metres 


*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  read¬ 
ings,  use  the  following  formula:  C  *  (5/9) (F  -  32).  To  obtain 
Kelvins  (K)  readings,  use:  K  -  (5/9) (F  -  32)  +  273.15. 
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AIR  DEMAND  TESTS,  LIBBY  DAM,  KOOTENAI  RIVER,  MONTANA 


PART  I:  INTRODUCTION 


Pertinent  Features  of  the  Project 

1.  Libby  Dam  (Figure  1)  is  located  on  the  Kootenai  River  in  north¬ 
western  Montana  221.9*  river  miles**  upstream  from  its  confluence  with 
the  Columbia  River  and  approximately  17  river  miles  above  the  town  of 
Libby,  Montana  (Figure  2).  It  is  a  multipurpose  project  constructed  as 


Figure  2.  Vicinity  map 


*  Based  on  official  river  miles.  Previous  publications  list  the  dam 
location  as  river  mile  219.0 

**  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (SI)  units  is  given  on  page  3. 
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an  integral  unit  of  the  comprehensive  water  resource  development  plan  of 
the  Columbia  River  Basin  in  the  United  States  and  Canada. 

2.  The  completed  project  consists  of  a  concrete  gravity  struc¬ 
ture,  a  reservoir  (Lake  Koocanusa)  having  a  total  gross  capacity  of 
5,870,000  acre-ft,  at  maximum  lake  elevation  2459  ft  msl,  and  a  hydro- 
power  installation  of  four  105,000-kw  units.  Four  additional  units  are 
now  under  construction.  The  maximum  length  of  the  lake  is  90  miles 
(42  of  which  extend  into  Canada).  The  lake  provides  4,980,000  acre-ft 
of  flood  storage. 


Outlet  Works 


3.  Desired  flow  through  the  structure  is  accomplished  by  means  of 
a  two-bay  ogee  spillway  over  which  flow  is  controlled  by  two  48-ft-wide 
by  56-ft-high  tainter  gates  and  three  sluices,  each  10  ft  wide  by  22  ft 
high,  controlled  by  10-ft-wide  by  17-ft-high  tainter  gates.  Both  the 
spillway  and  sluices  empty  into  the  same  hydraulic  jump-type  stilling 
basin  with  sloping  end  sill  and  no  baffles. 

Purpose  and  Scope  of  Tests 


Background 

4.  Libby  Dam  became  operational  in  March  1972.  Major  cavitation 
damage  occurred  in  the  center  and  right  sluices  in  September  1973  and 
July  1974,  respectively.  Following  repairs  to  the  center  sluice,  the 
U.  S.  Army  Engineer  District,  Seattle,  requested  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES)  Hydraulics  Laboratory  to  conduct 
field  measurements  to  determine  the  pressures  acting  on  the  invert  of 
the  center  sluice.  Tests  were  conducted  in  1974  and  a  report  of  find¬ 
ings  was  published  in  1976  (Hart  and  Tool  1976).  With  this  information 
as  reference  material,  WES  conducted  a  model  study  of  the  sluice  (Dortch 
1976).  As  a  result,  an  aerator  device  was  designed  that  would  aerate 
the  flow  along  the  sluice  boundaries  without  adversely  affecting  flow 
conditions . 
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5.  The  recommended  aeration  system  installation,  completed  by  the 
Seattle  District  in  the  summer  of  1980,  consists  of  a  slot  in  both  side- 
walls  and  the  floor  of  each  sluice  as  shown  in  Plate  1.  Observation 
during  initial  sluice  operation  revealed,  as  expected,  a  significant  in¬ 
crease  in  air  demand  above  the  premodified  condition.  The  Seattle  Dis¬ 
trict  requested  a  proposal  to  measure  the  air  demand  during  operation 

of  the  modified  sluices.* 

Purpose 

6.  The  center  sluice  of  the  project  receives  air  indirectly  from 
the  outer  sluice  air  vents  through  crossover  passages  (Plate  1) .  Tests 
were  conducted  to  determine  the  air  drawn  into  one  of  the  outer  sluices 
as  a  result  of  the  addition  of  the  aerator.  This  information  would  be 
used  to  design  an  independent  air  vent  for  the  center  sluice.  A  por¬ 
tion  of  one  of  the  aerators  is  shown  in  Figure  3. 

Scope 

7.  The  test  program  measured  air  demand  in  the  right  sluice  with 

a  full  range  of  gate  openings  at  1-ft  increments  from  1  through  16  ft 
with  essentially  full  pool  (2447.7  ft  msl) .  The  test  sluice  and  air 
vent  are  shown  in  Figure  4.  The  4-  by  4-ft  crossover  vent  to  the  center 
sluice,  shown  in  Plate  1  and  Figure  5,  was  blocked  to  prevent  airflow  from 
the  other  two  sluices.  Air-demand  measurements  included  the  air  sup¬ 
plied:  (a)  by  the  7-  by  5-ft  vent,  (b)  by  the  3.5-ft-diam  equipment 

shaft,  and  (c)  along  the  roof line  of  the  sluice  itself.  Separate  tests 
were  accomplished  with  the  equipment  shaft  open  and  closed.  In  addition, 
pressure  drops  from  the  atmosphere  to  the  sluice  chamber  were  measured; 
measurement  points  are  shown  in  Table  1  and  Plate  1. 


*  NPSEN-HH-HC,  letter  to  WES,  subject:  "Libby  Dam,  Prototype  Air- 
Demand  Tests,"  26  June  1980. 
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Figure  5.  Air  vent  system,  looking  upstream. 
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PART  II:  TEST  FACILITIES,  EQUIPMENT,  AND  PROCEDURES 


Test  Facilities 


Air  velocity 

8.  Pitot  tube  differential  pressures  were  measured  at  the  loca¬ 
tions  listed  in  Table  1  and  shown  in  Plate  1  for  determining  the  air 
velocities  in  the  air  vent,  equipment  shaft,  and  sluice.  For  calculating 
average  velocities,  the  average  differential  pressure  recordings  from  the 
oscillograms  were  used.  The  largest  pitot  tube  differential  pressure  ob¬ 
served  on  either  the  digital  multimeter  display  or  the  oscillogram  was  used 
in  calculating  the  maximum  recorded  velocity  at  each  pitot  tube  for  each 
gate  opening.  Typical  oscillograms  are  shown  in  Plates  2  and  3. 

9.  Pressure  probes,  such  as  the  pitot  tube,  rely  upon  a  unique 
relationship  that  exists  between  pressure  and  velocity.  This  relation¬ 
ship,  discussed  in  detail  in  cited  references  (Rouse  1962,  Bryer  and 
Pankhurst  1971),  is  presented  here  in  mathematical  form: 


where 

Ap  *  pressure  difference  between  positions  A  and  B 
pg  =  stagnation  pressure,  position  A 
Pq  *  pressure  of  the  undisturbed  flow,  position  B 
p  *  mass  density  of  air  at  ambient  temperature 
V  *  point  velocity  of  the  air  being  measured 

10.  Because  air  is  compressible,  higher  densities  and  tempera¬ 
tures  will  occur  at  the  pitot  tube  position  A,  the  point  of  maximum  pres¬ 
sure.  The  resulting  effect  on  the  differential  pressure  being  recorded 
is  directly  proportional  to  the  air  velocity.  For  Mach  numbers  less 
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than  unity,  this  effect  is  considered  negligible  (Vennard  1954).  The 
Mach  number  is  a  ratio  of  the  air  velocity  to  the  propagation  velocity 
of  a  pressure  wave  in  the  same  medium.  For  a  point  velocity  of  156  fps 
(see  Table  2,  Test  9),  the  Libby  tests  Mach  number  is  less  than  0.2. 
Therefore,  the  compressibility  effect  was  not  considered  in  converting 
the  measured  differential  pressures  to  point  velocites. 

11.  A  typical  arrangement  of  the  pitot  tube,  differential  pres¬ 
sure  transducer,  and  electrical  cables  is  shown  in  Figure  6.  The  trans¬ 
ducer  is  a  Validyne  DP9  differential  pressure  transducer  with  a  range  of 
0  to  0.5  psid.  Specifications  for  all  transducers  used  in  the  test  pro¬ 
gram  are  presented  in  Table  1.  Pitot  tubes  used  in  the  tests  were  cali¬ 
brated  in  a  wind  tunnel  at  the  National  Space  Technology  Laboratories 
(NSTL),  Bay  St.  Louis,  Mississippi.  The  calibration  provided  a  correc¬ 
tion  coefficient  for  the  pitot  tubes  used  in  the  study. 

12.  The  theoretical  coefficient  (k)  of  Equation  1  is  /2/p  .  The 


Figure  6.  Equipment  shaft  strut  (panel  removed)  and  instrumentation 
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average  temperature  during  the  field  tests  was  41°F  giving  a  value  of 
2  4 

p  =  0.00247  lb-sec  /ft  .  Therefore 


X  144 


00247 


=  341.47 


in. -ft 


sec 


/lb 


where  Ap  (Equation  1)  is  in  psi.  The  calibrated  pitot  tube  coefficient 
was  calculated  to  be  351.90  based  on  information  provided  by  NSTL.  A 
graph  comparing  the  coefficients  is  presented  in  Plate  4.  The  test  sec¬ 
tion  of  the  wind  tunnel  is  shown  in  Figure  7. 


i 


Figure  7.  Wind  tunnel  test  section 
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Struts 


13.  The  streamlined  struts  were  fabricated  from  mahogany.  A  typi¬ 
cal  strut  cross  section  is  shown  in  Figure  8.  The  cross  section  of  each 
consisted  of  two  ellipses  joined  at  their  minor  axes.  The  downstream 
major  axis  (a^)  was  longer  than  the  upstream  (a^) .  The  minor  axes  (b) 
were  identical.  The  struts  consisted  of  two  halves,  divided  along  the 
major  axis.  Sections  of  each  inner  side  of  the  strut  were  grooved  to 
accommodate  the  pitot  tubes,  pressure  transducers,  accelerometers,  and 
cables,  as  in  Figure  6. 


Strut  Location 

Cross 

-Section  Dimensions,  in 

al 

a2 

b 

Air  vent 

5.00 

8.50 

2.38 

Equipment  shaft 

4.00 

5.50 

1.75 

Sluice 

3.00 

3.00 

1.25 

Figure  8. 

Strut  cross 

section 
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14.  Strut  brackets  were  fabricated  to  attach  to  the  walls  and 


encase  the  strut  ends.  After  bolting  the  brackets  in  place,  shims  and 
setscrews  were  used  to  hold  the  strut  securely  in  place.  Electrical 
cables  that  exited  from  one  end  of  the  strut  were  attached  to  anchor 
bolts  every  foot  through  the  area  of  high  velocity.  Figure  9  shows  the 
air  vent  strut  in  place  with  brackets  and  secured  cables. 


Figure  9.  Air  vent  strut  and  pitot  tubes 


Transducer  locations 

15.  The  pitot  tube  and  acceleromenter  arrangements  are  presented 
in  Table  1  and  Plate  5.  The  equipment  shaft  pressure  transducer  (ESI) 
and  accelerometer  (ESA)  were  located  in  the  center  of  the  strut.  The 
distance  ratios  (y/r)  from  the  air  vent  strut  center  line  to  the  pitot 
tubes  are  shown  in  Plate  5.  The  sluice  strut  was  attached  to  the  sluice 
roof  as  shown  in  Figure  10  and  Plate  1.  The  pitot  tubes  were  located 
6  and  13.5  in.  from  the  roof;  this  was  necessary  to  keep  them  dry  during 
the  higher  gate  openings. 


Figure  10.  Sluice  strut  and  pitot  tubes 
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Strut  acceleration 

16.  Computations  indicated  the  possibility  of  the  struts  vibrat¬ 
ing  at  or  near  their  natural  frequency.  For  this  reason,  a  +5  g  accel¬ 
erometer  was  placed  in  the  center  of  the  air  vent  and  equipment  shaft 
struts  (AVA  and  ESA,  respectively) .  The  accelerometers  were  secured 
with  wood  screws  and  their  cables  led  out  of  the  struts  through  the  pres¬ 
sure  transducer  cable  passageway  grooves. 

Differential  pressure 

17.  Pressure  in  the  gate  ?b»sher  was  measured  with  a  differential 
pressure  transducer  (GC1,  Tsblo  j  .  ''  ..ate  1).  The  purpose  was  to  measure 
the  drop  in  pressure  from  4J  .*  vphere  to  the  aerator  which  caused  the 
airflow  through  the  air  venr*  lament  shaft,  and  sluice.  The  trans¬ 
ducer  was  placed  on  the  haiiv.  >  of  the  gate  chamber  walkway  a  short 
distance  from  the  aerator. 


Other  Measurements 


18.  Other  recorded  data  consisted  of  upper  pool  elevation,  air 
temperature,  barometric  pressure,  gate  opening,  and  water  discharge. 

These  data  were  provided  by  project  and  District  personnel.  Water  dis¬ 
charge  was  determined  from  computed  discharge  rating  curves. 

Recording  Equipment 

19.  The  recording  equipment  consisted  of:  (a)  WES-fabricated 
model  03  bridge  amplifiers  to  condition  the  strain  gage  transducers 
(GC1,  AVA,  ESA)  output  signal,  (b)  a  WES-fabricated  channel  selector, 

(c)  a  Validyne  CD12  recorder  for  providing  an  analog  output  signal, 

(d)  a  Fluke  8300A  digital  multimeter  for  displaying  the  Validyne  re¬ 
corder  output  voltage  digitally  (AVI,  AV2,  AV3,  AVA,  ESI,  SL1,  and  SL2) , 
and  (e)  a  CEC  model  5-12A,  7-in.  direct  print  oscillograph  capable  of 
reproducing  up  to  18  channels  of  data  at  chart  speeds  of  1/A  to  6A  ips. 
Chart  speeds  of  A,  16,  and  6A  ips  wer?  used  during  the  Libby  tests.  Fig¬ 
ure  11  shows  the  equipment  setup  at  the  recording  station. 
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Figure  11.  Test  recording  equipment 
Test  Procedures 

20.  Tests  were  conducted  on  22  October  1980.  Data  for  each  gate 
opening  were  recorded  on  the  oscillograph  chart  and  also  read  from  the 
multimeter  and  recorded  manually.  The  procedure  was  the  same  for  all 
tests  and  consisted  of  the  following: 

_a.  Record  step  calibrations. 

_b.  Record  test  number,  date,  time,  and  gate  opening. 

£.  Raise  gate  to  desired  opening;  allow  flow  to  stabilize. 
d_.  Set  channel  selector. 

e.  Set  zero  and  gain  on  Validyne  recorder. 

f^.  Record  data  on  the  oscillograph  chart  and  manually. 

Repeat  steps  d-f  until  all  data  channels  are  recorded. 

h.  Record  upper  pool  elevation,  barometric  pressure,  and 
air  temperature. 

_i.  Record  step  calibrations. 


Notes  were  made  on  the  oscillograms  during  the  tests  for  reference 
during  data  reduction  and  analysis. 


PART  III:  TEST  RESULTS  AND  ANALYSIS 

Air  Discharge 

21.  Point  velocities  in  the  equipment  shaft,  air  vent,  and  sluice 
were  calculated  by  the  methods  previously  described.  Conditions  at  each 
measurement  station  were  different,  necessitating  different  methods  of 
coverting  the  calculated  point  velocities  to  average  velocities  in  their 
respective  conveyance  system.  These  methods  and  results  are  described 
below.  In  addition,  the  pressure  drop  data  measured  at  the  gate  chamber 
were  used  to  approximate  losses  through  the  equipment  shaft  and  air  vent 
Equipment  shaft 

22.  Rouse  (1962)  shows  that  when  the  Karman-Prandtl  velocity  dis¬ 
tribution  equations  for  smooth  and  rough  surfaces  are  each  solved  for 
the  mean  velocity  and  then  subtracted  from  their  original  form,  an  iden¬ 
tical  expression  results.  With  a  slight  adjustment  of  coefficients,  the 
derived  equation  was  found  to  agree  very  closely  with  experimental  data. 
This  final  form,  which  presents  the  mean  velocity  (V)  as  a  function  of 

the  maximum  (center  line)  velocity  (V  )  and  the  resistance  coefficient 

m 

(f),  is 

V 

V  - - - -  (2) 

1.43  /f  +  1 

23.  Provided  the  Reynolds  Number  (St)  is  sufficiently  high,  the 
velocity  distribution  in  a  circular  duct  should  be  typical  of  fully  de¬ 
veloped  turbulent  flow  at  a  distance  of  about  50  diameters  or  more  down¬ 
stream  of  the  entrance.  The  equipment  shaft  strut  was  located  49  diam¬ 
eters  downstream  of  the  entrance  (top  of  the  dam)  and  SI  was  equal  to  or 
greater  than  7  x  10^  (Tables  2  and  3).  Therefore,  Equation  2  is  appli¬ 
cable  for  converting  the  measured  center-line  velocities  to  mean  equip¬ 
ment  shaft  velocities. 

24.  The  Corps  of  Engineers  (OCE  1980)  recommends  a  surface  rough¬ 
ness  height  (k)  value  of  0.0020  for  circular  concrete  conduits.  For  a 
3.5-ft-diam  conduit,  this  corresponds  to  a  resistance  coefficient  (f) 
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value  of  approximately  0.019.  Since  this  is  a  design  recommendation,  a 
further  review  was  conducted.  The  k  value  range  for  concrete  is  given 
in  most  references  as  0.001  to  0.01  ft.  This  corresponds  to  f  values 
for  this  study  of  0.015  and  0.026,  respectively.  Albertson  et  al.  (1961) 
state  "the  average  value  of  k  should  be  utilized  unless  additional  in¬ 
formation  is  available."  The  average  value  of  f  then  is  0.0205. 

25.  A  value  of  f  =  0.021  was  used  in  the  computations.  It  is 

noted  that  using  the  above  extreme  values  of  f  and  a  center-line 

velocity  of  100  fps,  the  average  velocity  differs  by  only  A  percent. 

The  flow  rate  was  then  determined  by  multiplying  the  average  velocities 

2 

by  the  equipment  shaft  area  (9.62  ft  ).  Plate  6  presents  the  average 
equipment  shaft  discharge  values  for  each  gate  opening.  Table  A  also 
presents  the  average  as  well  as  maximum  computed  discharges.  Tables  2 
and  3  present  the  observed  maximum  and  average  equipment  shaft  point 
velocities  which  were  used  in  the  air  discharge  computations. 

Air  vent 

26.  Point  velocities  were  measured  at  four  vertical  locations  in 
the  air  vent  as  shown  in  Plate  5.  The  strut  was  located  57  ft  down¬ 
stream  of  the  vent  entrance  or  9.8  equivalent  diameters  (D^  =  5.75  ft). 
The  strut  was  situated  in  this  forward  position  in  an  attempt  to  mini¬ 
mize  back-pressure  effects  in  the  relatively  short  vent  (15. A  D^)  which 
exits  at  a  right  angle  into  the  gate  chamber. 

27.  Because  of  the  short  distance  from  the  entrance  (and  after 
reviewing  the  data),  the  velocity  distribution  was  assumed  to  be  essen¬ 
tially  constant  from  wall  to  wall.  The  mean  velocity,  then,  was  assumed 
to  be  the  average  value  of  the  four  measurements.  To  determine  the 
degree  of  validity  of  this  assumption,  the  standard  deviation  (o)  of 

the  velocities  for  each  gate  opening  was  computed.  The  average  of  this 
value  divided  by  the  average  of  the  means  (V)  gives  the  coefficient  of 
variation  (CV)  (Miller  and  Freund  1977),  i.e., 

CV  =  -  (3) 

V 

28.  For  the  Libby  air  vent  measurements,  the  value  of  CV  was 
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determined  to  be  0.080  and  0.095  for  the  average  and  maximum  recorded 
point  velocities,  respectively.  In  other  words,  on  the  average,  the  re¬ 
spective  standard  deviations  were  8.0  and  9.5  percent  of  the  mean.  This 
implies  that  the  assumption  of  uniform  velocity  distribution  at  the  strut 

is  reasonable.  These  mean  velocities  were  multiplied  by  the  air  vent  area 
2 

|  (34.5  ft  )  to  determine  the  discharge.  Air  vent  discharge  versus  gate 

opening  is  presented  in  Plate  6  with  the  equipment  shaft  open  and  closed. 
Note  that  the  discharge  is  approximately  the  same  in  each  case.  Table  4 
presents  the  average  and  maximum  computed  air  vent  discharges.  The  maxi¬ 
mum  values  are  approximately  20  percent  and  25  percent  greater  than  the 
average  values  with  the  equipment  shaft  open  and  closed,  respectively. 
Tables  2  and  3  present  the  observed  maximum  and  average  air  vent  point 
velocities  which  were  used  in  the  air  discharge  computations.  The  four 
poitit  velocities  used  in  each  computation  were  not  recorded  simultaneously. 
Plate  7  presents  the  air  to  water  discharge  ratio  as  a  function  of  the 
Froude  number  plotted  on  the  Hydraulic  Design  Criteria  Chart  050-1. 

29.  Sudden  gusts  of  high-velocity  air  occurred  during  many  gate 
openings.  To  estimate  the  relative  difference  between  a  gust  velocity 
and  the  recorded  high  velocities  the  gust  peak  of  AV4,  test  26  (Table  2 
and  Plate  2),  was  compared  with  the  average  of  the  recorded  maximum  velo¬ 
cities  of  AVI,  AV2,  and  AV3,  same  test.  These  values  are  139.0  and  103.0 
fps,  respectively.  This  indicates  that  the  gust  velocity  was  35  percent 
greater  than  the  average  of  the  recorded  nongusting  maximum  velocities. 

S3  uice 

30.  After  the  aerators  were  installed,  spray  was  observed  being 
drfcwn  back  into  the  sluices  at  their  exit  Dortals.  Therefore  it  was 
determined  that  airflow  in  the  sluice  should  be  monitored  also.  The 
sljuice  air  velocity  was  measured  near  the  roof  as  shown  in  Plate  1.  As 
stated  previously,  this  was  necessary  to  ensure  that  the  strut  and  encased 
transducers  cleared  the  discharging  water  at  all  gate  openings.  It  was 
assumed  that  a  bidirectional  velocity  profile  existed  in  the  sluice,  i.e., 
the  pressure  differential  causing  air  near  the  roof  to  move  upstream  and 
shear  at  the  air-water  interface  causing  downstream  flow  of  the  lower 
portion.  This  is  similar  to  the  theoretical  study  of  flow  between  two 
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parallel  plates  (Schlichting  1968).  In  this  case,  one  wall  is  at  rest 
(sluice  roof)  and  the  other  wall  moves  at  velocity  U  (water-surface 
velocity) . 

31.  Use  of  this  theoretical  approach  for  computing  sluice  dis¬ 
charge  was  complicated  by  the  fact  that  the  velocity  data  (SL1,  SL2) 
indicated  that  (a)  there  was  very  little  airflow  at  the  measurement 
point  for  gate  openings  of  11  ft  and  less,  and  (b)  the  measured  flow 
oscillated  in  both  directions  at  the  higher  gate  openings.  Figure  12 
presents  an  illustration  of  this  phenomenon.  For  these  reasons,  it  is 
believed  that  very  little  of  the  air  entering  the  sluice  exit  actually 
reached  the  aerator,  instead  was  drawn  back  out  by  shear  along  the  air- 
water  interface. 
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TEST  22.  GATE  OPENING  11  FT 
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TEST  20.  GATE  OPENING  13  FT 


Figure  12.  Sluice  airflow  (SL2) 


Pressure  Differential 


32.  The  pressure  differential  between  the  atmosphere  and  the  gate 
chamber  was  measured  with  a  differential  pressure  transducer  located 
near  the  aerator  (GC1),  which  is  shown  in  Plate  1.  These  data  were  used, 
with  the  velocity  measurements,  to  approximate  losses  through  the  air 
vent  and  equipment  shaft.  The  Bernoulli  equation  (Rouse  1962)  may  be 
used  to  evaluate  the  total  energy  loss  from  the  entrances  to  the  gate 
chamber  as  follows: 

2  2  2 
V  V  V 

P  r  +  Pe  *  P  /  +  Pc  +  P  f  K 


where 


p  *  density  of  air,  slugs/ft' 
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V  =  entrance  velocity,  ft/sec 

e  2 
Pg  =  entrance  pressure,  lb/ft 

V  =  chamber  velocity,  ft/ sec 

C  2 
Pc  =  gate  chamber  pressure,  lb/ft 

V  =  measurement  point  (air  vent  or  equipment  shaft)  velocity, 
m  ft/sec 

K  =  loss  coefficient 

Assuming  Vg  <<  <<  Equation  4  reduces  to  Ap  =  p 

2A  ^ 

°r  K'  ^ 


where 


K'  =  A 

m  \  pK 


where 


A^  =  conveyance  system  cross  section  at  the  measurement  location 

Ap  *  pressure  differential  between  the  atmosphere  and  the  gate 
chamber 

33.  Using  the  field  data,  the  values  of  ,  v^Ap  ,  and  K' 
(Equation  5)  for  each  gate  opening  were  computed  for  both  the  equipment 
shaft  and  the  air  vent  flows.  A  plot  of  the  individual  data  points  with 
their  corresponding  averaged  K’  value  is  presented  in  Plate  8.  With 
this  information  the  apparent  total  loss  coefficient  (K)  through  each 
system  can  be  determined  with  Equation  6.  For  the  equipment  shaft  and 
air  vent,  the  values  are: 


2(9.62) 


(0.00247)  (176.0)' 


2(34.5) 

(0.00247)  (572. 1): 


=  2.9 


34.  These  K  values  are  the  sum  of  the  loss  coefficients  for 


each  conveyance  system.  The  primary  individual  losses  contributing  to 
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the  total  were  assumed  to  be  the  entrance,  friction,  and  exit  losses. 
Assigning  values  of  0.5,  0.02  L/D  (conveyance  system  length  to  diameter 
ratio),  and  1.0  to  these  losses,  respectively,  the  sums  were  determined 
to  be: 


Kes  -  2.59  Kav  -  1.82 

The  equipment  shaft  coefficient  seems  reasonable  while  the  air  vent 
value  is  lower  than  the  computation  based  on  the  test  data.  The  addi¬ 
tional  losses  may  be  due  to  the  combined  sudden  expansion,  90-deg  turn, 
and  grating  (el  2265.63,  Plate  1)  at  the  top  of  the  gate  chamber. 

Strut  Vibrations 


35.  The  aformentioned  accelerometers  were  placed  in  the  equipment 
shaft  and  air  vent  struts  to  monitor  their  vibrations.  Prior  to  the 
tests,  each  was  struck  lightly  with  the  butt  of  the  hand  and  the  re¬ 
sponse  recorded.  These  recordings  are  shown  in  Figure  13.  During  the 
entire  test  program,  the  air  vent  strut  (AVA)  vibrated  very  near  its 
natural  frequency.  The  maximum  recorded  peak-to-peak  displacement  was 

_3 

5  x  10  in.  at  a  frequency  of  56  Hz  and  a  gate  opening  of  10  ft.  The 

-4 

equipment  shaft  strut  (ESA)  maximum  displacement  was  7.5  x  10  in.  at 
130  Hz  (9-ft  gate  opening). 


ESA  fn  =  183  Hz 


Figure  13.  Air  vent  strut  (top)  and  equipment  shaft 
strut  (bottom)  "tap"  response 
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Measurement  Errors 


36.  The  maximum  error  in  pressure  differential  determination  from 
the  records  or  the  digital  multimeter  is  estimated  to  be  approximately 
7  percent.  This  value  is  based  on  the  following  estimated  errors. 

a.  The  step  calibrations  were  read  to  within  +0.02  in.  over 
a  1-in.  deflection  for  an  error  of  2  percent. 

t>.  The  oscillogram  traces  were  read  to  within  +0.02  in.  over 
a  deflection  of  0.4  in.  for  an  error  of  5  percent. 

£.  It  is  estimated  that  the  digital  multimeter  was  read  with 
an  error  of  +5  percent. 
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PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


37.  The  following  determinations,  conclusions,  and  recommenda¬ 
tions  result  from  literature  review,  field  observations,  and  analysis  of 
the  reduced  data  of  the  Libby  Dam  air  demand  tests. 

a..  Measurements  indicate  that  practically  all  of  the  air 
is  drawn  into  the  aerator  through  the  air  vent  and 
equipment  shaft  (when  open). 

Ij.  The  air  vent  draws  approximately  the  same  amount  of  air 
whether  the  equipment  shaft  is  open  or  closed. 

c^.  The  maximum  measured  airflow  values  of  the  air  vent  ex¬ 
ceeded  the  average  values  by  approximately  20  and  25  per¬ 
cent  with  the  equipment  shaft  open  and  closed, 
respectively. 

cL  Observations  were  made  of  air  being  drawn  into  the  exit 
portal  of  the  sluice  at  many  gate  openings.  Air  in  the 
sluice  oscillated  in  both  directions  at  the  measurement 
point,  indicating  it  may  not  have  reached  the  aerator. 

A  reasonable  approximation  of  the  loss  coefficients 
through  the  equipment  shaft  and  air  vent  is  2.5  and  2.9, 
respectively. 

_f.  Since  the  spillway  location  precludes  duplicating  the 
existing  air  vent  geometry,  a  loss  coefficient  for  the 
proposed  center  sluice  air  vent  cannot  be  approximated 
directly  from  these  measurements.  If  the  surface  finish 
is  similar  to  the  equipment  shaft  and  air  vent  tested, 
0.02  would  be  a  reasonable  approximation  of  the  friction 
coefficient  to  be  used  in  estimating  friction  losses, 

j>.  At  most  gate  openings  the  airflow  rate  was  not  constant, 
experiencing  sudden  gusts  of  high  velocity.  This  insta¬ 
bility  should  be  considered  in  the  air  vent  design. 

h.  Because  of  j>,  future  air  demand  tests  should  consist  of 
simultaneous  recording  of  all  channels  on  magnetic  tape 
in  order  to  fully  analyze  the  gusting  phenomenon. 

jl.  Model  results  (Dortch  1976)  indicate  that  an  increase  in 
air  demand  will  occur  during  combined  spillway-sluice 
flow.  This  should  be  considered  if  combined  flow  is 
anticipated. 
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Scaled  from  drawings. 


Air  Vent  and  Equipment  Shaft  Maximum  Recorded  Point  Velocities 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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Cover  Title. 
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